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Abstract 



To explain the experimental spectra for angle resolved photoemission we con- 
J^ i sider a modified t — J model. The modified model includes next nearest (t') 

and next next nearest (t") hopping as well as Hubbard model corrections to 
the spectral weights. A Dyson equation which relates the single hole Green's 



^T^ . functions for a given pseudospin and given spin is derived and is compared to 

-y ' experimental results. 
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I. INTRODUCTION 

Recent angle resolved photoemission (ARPES) measurements by Wells et aB and by 
Pothuizen et aB for insulating Copper Oxide Sr2Cu02Cl2 provide a unique possibility to 
experimentally determine the single hole properties. In the frameworks of a t — t' — J model 
the problem has been analyzed by Nazarenko et aB using a cluster method and Bala, Oles, 
and ZaanenQ using a self-consistent Born approximation (SCBA). It was demonstrated in 
these works that including t' hopping improves agreement with experiment, but a complete 
description of the ARPES spectra especially the k dpendent intensity was not achieved. 

In the present work we consider a, t — t' — t" — J model. A Dyson equation which relates 
the one electron Green's function measured in experiment to the hole Green's function found 
in the self-consistent Born approximation is derived. We also introduce into this equation 
the corrections which originate from finite U as in the Hubbard model. With parameters 
taken from LDA band calculations by Andersen et a£l we got a reasonably good description 
of the experimental ARPES spectra . The results are sensitive to t" . The parameter t' is 
less important. The importance of t" was also pointed out by Belinicher and ChernyshevQ. 
For the explanation of width of ARPES spectra we need to introduce additional broadening 
and possible origins of this are discussed . 

The Hamiltonian of t — t' — t" — J model is of the form 

H = —t 2_^ c^^Cju — t 2_^ Cj^Cj^o- — t 2^ c^^Cj^a + "^ Z^ SjSj. (1) 

(ij)a {ijl)o- {m)a {ij)cr 

c\„ is the creation operator of an electron with spin a (a =t, i) at site i of the two-dimensional 
square lattice, the {ij) represents nearest neighbor sites, {iji) - next nearest neighbor (diago- 
nal), and {iJ2) represents next next nearest sites. The spin operator is Sj = |c|Q,o"a/3Cj^. The 
size of the exchange measured in two magnon Raman scatteringtMa is J = 125meV. The 
most recent calculation of the hopping matrix elements has been done by Andersen et aE. 
They consider a two-plane situation and the effective matrix elements are slightly different 
for symmetric and antisymmetric combinations of orbitals between planes. After averaging 
over these combinations we get: t = 386meV, t' = —105meV , t" = 86meV. Below we set 
J = 1, in these units 

t = 3.1, t' = -0.8, t" = 0.7 (2) 



II. HOLE GREEN'S FUNCTION WITH FIXED PSEUDOSPIN. 
SELF-CONSISTENT BORN APPROXIMATION (SCBA) 

It is well known that at half filling (one electron per site) the model under considera- 
tion is equivalent to a Heisenberg model. It represents a Mott insulator with long range 
antiferromagnetic order. We denote the corresponding ground state wave function by |0). 
We are interested in the situation when one electron is removed from this state, so a single 
hole is produced. The dynamics of a single hole in the antiferromagnetic background can be 
described by SCBABI12I. This approximation works very well due to the absence of a single 
loop correction to the hole-spin-wave vertexlii]tJ. Let us recall the idea of this approxima- 
tion. The bare hole operator di is defined so that d] oc Cj| on the | sublattice and oc Cj| on 
the I sublattice. In the momentum representation 



4 = /^ Ya.e^^^^, dl^ = .—-^ yc.^e"^^^. (3) 

N is number of sites, m = KOl^i^lO)! ~ 0.3 is the average sublattice magnetization . The 
quasi- momentum k is limited to be inside the magnetic Brillouin zone: 7k = ^{coskx + 
cos ky) > 0. In this notations it looks like d^^ has spin a = ±1/2, but actually rotation 
invariance is violated and a is a pseudospin which denotes the sublattice. Nevertheless the 
pseudospin gives the correct value of the spin 2;-projection: Sz = <J = ±1/2. The coefficients 
in (|^) provide the correct normalization: 

The hole Green's function is defined as 

Gd{e,k) = -I j {Q\Td^,{T)diM\^)d^^dT (5) 

The t', t" terms in the Hamiltonian (|l]) correspond to the hole hopping inside one sublattice. 
This gives the bare hole dispersion 

eok = 4t' cos kx cos ky + 2t"(cos 2kr^ + cos 2ky) -^ /9oi7k ± /3o2(7k )^, (6) 

where 7,^ = \{cos k^ — cos ky), /?oi = 4(2t" + t'), and /?02 = 4(2t" — t'). In equation (H) we took 
into account that the sign of a hole dispersion is opposite to that for an electron (maximum 
of electron band correspond to minimum of hole band), and omitted some constant. The 
bare hole Green's function is 



Go,(e,k) = — -. (7) 

e - eok + ^0 

For spin excitations the usual linear spin-wave theory is used (see, e.g. review papertj). 

It is convenient to have two types of spin- waves, aj, with S^ = —1, and M with Sz = +1, 

and q restricted to be inside the magnetic Brillouin zone. 



^^'- 



e-''!"' ^ UqCtq + t;q/5Lq, (8) 






The spin-wave dispersion and parameters of Bogoliubov transformation diagonalizing the 
spin- wave Hamiltonian are: 






^^q = -signi'y^) 



CUq 2' 



Hopping to nearest a neighbor in the Hamiltonian ([l|) gives an interaction of the hole 
with spin-waves. 

Hh,sw = Yl 9KQ (4+qiC?kTaq + 4+qT^ki/5q + H.C.^ , (lO) 

k,q 

with vertex 5fk,q given by 

9k,q = (0|aqC?kT|i^t|4+qJO) = 

2 ^ e-'^'-^+^(''+'^>»((0|aqc]^c,TC,V,T|0) + (0|aqc]^c,icl^Qt|0) 



iV(l/2 + m) ,-^^^ 



2 



E e-'''-^+^(''+^)^^(0|aq (57 + Sr) |0) = ^t^^il^u^ + Tk+qt^q). (11) 



In this calculation we have used the usual mean field ground state factorization approxi- 
mation: (0|aqc]|Cj|q|Cj^|0) ~ (0|aqc]|Cj||0)(0|cJ|Cj||0) = (0|Q;q5'7|0)(l/2 + m). The vertex 
(7k,q is independent of t', t" because these parameters correspond to hopping inside one sub- 
lattice. The form of the vertex (7k,q is well known. The actual purpose of calculation (pUf ) 



is to demonstrate that (7k,q is valid in a more general situation than it is usually believed. 
It is independent of the particular value of the sublattice magnetization m. Therefore, for 
example, f7k,q remains the same in the presence of strong additional frustrations. 

One can easily prove that the spin structure of the interaction ( JIDD forbids single loop 
corrections to the hole-spin-wave vertex and, as usually, the two loop correction is small 
numericallyHj cia. So, due to the spin structure we have an analog of the well known Migdal 
theorem for electron-phonon interactions. This justifies SCBA according to which the hole 
Green's function satisfies a simple Dyson equation 

Gd{e, k) = f e - eok - E ^k-q,qG'd(e -uj^,k-q)+to] . (12) 

The anomalous Green's function —i{0\Tdki{t)dl^^{Q)\Q) vanishes because the z-projection of 
spin is conserved. Due to the definition of the operators (|^) the Green's function (|]) is 
invariant under translation with the inverse vector of the magnetic sublattice Q = (±7r, ivr) 

G,(e,k + Q) = G',(e,k). (13) 

The numerical solution of equation (|T2|) is straightforward. As usual, to avoid poles we 
replace zO — >■ ir/2 = i 0.1. The energy scale consists of 300 points with variable density 
(concentrated near sharp structures of Gd)- The number of points in the magnetic Brillouin 
zone is 10^ which is equivalent to the lattice 140 x 140. The plots of — - Im G'd(e, k) as a 
functions of e for k = (7r/2,7r/2), k = (7r/2,0), k = (7r,0), and k = (0,0) are presented in 
Fig.l. We recall that we use the set of parameters (0) based on Ref.H. The position of the 
lowest peak gives the quasiparticle energy. Results of the calculation can be fitted by the 
formula 

ek = const + (3,^1 + /?2(7k )' + ^2(1^)' (14) 

Pi ^ 3.0, P2 ~ 3.8, P'^ ^ -1.5. 

The dispersion has minima at (ko = (±7r/2, ±7r/2). The hole pockets are slightly stretched 
along the direction to the zone center, and it is very different from a pure t — J modeltlrE3. 
The quasiparticle residue Z^ can be found as the area under the peak. At the dispersion 
minimum it equals Z^ = 0.38. So it is larger than in a pure t — J modelcJ'E3, but away 
from the dispersion minimum it drops down very rapidly. Plots of the residue Z^ as a 
function of k for k G [(7r/2, 7r/2) - (0,0)] and k G [{71/2, 7t/2) - (0,7r)] are given at Fig. 
2. Throughout the Brillouin zone the residue is fitted by 

Zt = < (1 - 0.97^ - 1.52(7k )' + 0.057^ + 0.69(7k )' + 0.5(7k7k )') • (15) 



III. HOLE GREEN'S FUNCTION WITH FIXED SPIN. DYSON EQUATION 
RELATING TWO GREEN'S FUNCTIONS 

The operators c?kt, (ikj discussed in the previous section are defined at different sublat- 
tices. However, when a photon kicks out an electron from the system it does not separate 
the sublattices. Therefore for this process we have to define the operator as a simple Fourier 
transform 

Ck'^ = \/fEc-e^'"- (16) 

* i 

The normalization is chosen in such a way that (0|cj^|Ck||0) = -|:(0| J2i c1t^«t|0) ~ ^- ^^ ^^^ 
consider Ck^ as an external perturbation, and the corresponding Green's function is 

Ge(e,k) = -zy'(0|rcL(r)ck.(0)|0)e*^^rfr. (17) 

This is the Green's function measured in ARPES. Let us now find the relation between 
Gcle, k) and (^^(e, k). 

The operator Cko- acting on the vacuum (ground state of the Heisenberg model) can 
produce a single hole state. We denote the corresponding amplitude by Ok and show it in 
Fig. 3a as a cross. The thick line corresponds to the Green's function Gc (|T7|) and the thin 
line corresponds to the Gd (||). The amplitude Ok equals 

ak = (OMkTCkilO) = (0| [^ ^(i/2 + m) g^^V-"^^^) (/l?'^^'^'") '°^ = VV2 + m. 

(18) 

The operator Cko- acting on the vacuum can also produce a hole + spin-wave state. This 
amplitude is shown in Fig. 3b as a circled cross with the dashed line being a spin-wave. We 
denote this amplitude by 6k,q 



6.,. = (o|/?,4-.iCk,|o) = (o|/?, I J ^(i/2 + ^) g"'V'^^'"^"j f ylE^i^^"^^ I lo) 
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mJT.Ste'A\0) = ^v^. (19) 



We stress that (^) is a bare vertex. It corresponds to the instantanious creation of a hole + 
spin wave, but not the creation of a hole with a subsequent decay into hole + spin-wave. To 
elucidate this point look at Fig.4.The upper part of this figure describes the wave function 



of the initial Neel state: a - component without spin quantum fluctuations, b - component 
with spin quantum fluctuation. The lower part of Fig. 4 arises imediately from the upper 
one after kicking out an electron with spin up. Part a does not contain a spin flip, and it 
corresponds to the amplitude Ok (|TB]). Part b does contain a spin flip, and it corresponds to 
the amplitude 6k,q (|lD- Note that 6k,q ^ cxd at q ^ 0. The reason is that the operator (p!6D 
does not correspond to any quasiparticle of the system with long-range antiferromagnetic 
order, and therefore the usual Goldstone theorem is not applicable. 

Let us denote by a dot (Fig. 3c) the usual hole-spin- wave vertex f^k.q given by eq. (^T]). 
In the leading in t approximation the amplitude of single hole creation by the external 
perturbation ( [TB| ) is given by diagrams presented at Fig. 5, with the thin solid line in this 
case the bare hole Green's function (^. If we set t' = t" = and e = eok = calculation of 
this amplitude can be easily done analytically 

M,(e = 0, k) = ak + E ^'^""^-"^ = Vo:8 - ^ E '''^^'"'^'''^ ^ ^'"'^^ " 



q eOk — eok-q — ^q "^ q ^' 



q 



0-8 + ^7k E (^ - ^) = v^ (1 + 0-45 ■ t ■ 7k) . (20) 



N ^Xuj,, 2 

Mf is the quasiparticle residue of the Green's function ([171) . The eq. ( pOD agrees with the 
result obtained using a string represent at ionli^. We stress that even at t = the residue is 
0.8 due to the spin quantum fluctuation in the ground state of the Heisenberg modelE^O 

Now we can find the relation between Green's functions Gc ( p^ and Gd (|^)- In the 
leading in t approximation it is given by diagrams presented at Fig. 6 with the thin solid 
line being in this case the bare hole Green's function God (0)- ^ow let us dress these diagrams 
by higher orders in hopping t. As we already discussed there is no single loop correction to 
the "dot". We neglect double loop correction to the "dot" as it has been done in SCBA. 
Therefore the only possibility is an introduction of a self energy corrections. An example 
of the correction to diagram Fig. 6c is shown at Fig. 7. To take into account all these 
corrections we need just to replace at Fig. 6 all bare hole Green's functions (|^) by dressed 
hole Green's function given by eq. (jT^). So, the Fig. 6 actually represents a Dyson equation 
relating Gc (|17D and Gd (|5|). In analytical form it is 

Gc(e, k) = aiGd{e, k) + E b^c^Gdie - tUq, k - q) + 



+ 2akGrf(e, k) 



G',(e,k) 



E ^k,q5'k-q,qG'rf(e - CUq, k - q) 

L q 

n2 



+ (21) 



E ^k,qfi'k-q,qG'd(e - CUq, k - q) 



So as soon as we have found Gd using SCBA (jT^) we can calculate the Green's function Gc 
defined by eq. (|T^). The imaginary part of Gc{e, k) gives directly the spectra measured in 
ARPES experiments. 

We now discuss sum rules. All singularities of Green's functions are in the lower half 



plane of complex e. Therefore if we integrate eq.([T^) over e from — oo to +00, this integral 

can be replaced by the integral over an infinite semi-circle in the upper e half plane. For 

infinite e, Gd = God, and we get the well known sum rule 

1 r°° 
- -Im Gd(e,k)rfe = l, (22) 

TT J-00 

which agrees with with eq.@. If we integrate now eq. (|2lD in the same limits, the terms 
which contain more than one Green's function give zero contribution, because the integral 
can be transfered into the upper complex e half plane. And we find 

-^Jm|~ G,(6,k)rf6= (-i/m I G,(e,k)rfe) (^ + eO = 0-8 + | E< = 1- (23) 

Thus the equation (|2TD reproduces the correct normalization: (0|ci^|Ckt|0) = 1. 

The vertex 6^ q (0) is invariant under translation with the inverse vector of magnetic 
sublattice Q = (±7r, ±7r): &k+Q,q = ^k,q- At the same time the vertex ^f^q (|ll|) changes 
sign with this translation: (j'k+Q.q = — 5'k,q- Therefore the diagrams Fig. 6c, d change sign at 
k — i> k + Q and 

a(e,k+Q)7^G'e(e,k). (24) 

Due to the same properties of vertices 6k,q and f7k,q the diagrams Fig. 6c,d,e (square brackets 
in eq. (|2T|)) vanish at the face of magnetic Brillouin zone (7k = 0). The diagram presented 



at Fig. 6b (term with 6^ in eq. (^T]) ) is small numerically. Therefore at the face of magnetic 
Brillouin zone Gc(e, k) k, Gd{^, k). However away from the face they differ significantly. The 
plots of —^ Im Gc{e, k) as a functions of e for k = {n/2, vr/2), k = (7r/2, 0), k = (tt, 0), and 
k = (0, 0) are presented at Fig.8. A plot of the quasiparticle residue Z^ as a function of k 
along (1,1) direction is given at Fig. 9. The quasiparticle residue outside the magnetic zone 
is smaller than that inside. For comparison we also present a plot of Z^. 

IV. HUBBARD MODEL CORRECTION 

The picture considered above corresponded to a modified t — J model. It means that 
double electron occupancy was forbidden. Now we want to take into account the fact that 
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the t — t' — t" — J model originates from the Hubbard modeL We assume that it is a simple 
one band Hubbard model with on site repulsion U . First of all this gives some corrections to 
the "bare" hole dispersion (|^), see, e. g. Ref.u. However we assume that renormalization is 
done and these corrections are already included in the values of effective hopping amplitudes 
t', t" given in (|]). There are also some corrections to the hole-spin- wave vertexQ, but they 
are small at t', t" ^ U . The really important effect is the renormalization of the vertex Ok 
(p!8|). In t — J model this vertex is given by the process shown at Fig. 10a: an electron is 
removed from corresponding sublattice. In Hubbard model there is an additional possibility 
shown at Fig. 10b: first the electron hops to occupied nearest site, and then it is removed 
from this site. Simple calculation shows that this gives 

ak ^ ak X (^1 + -7k)) = v/l/2 + m (^1 + --f^ , (25) 

4t \ [Y / J 



&k,q ^ Ok,q X (^1 + — 7kj = ^—V^ (^1 + y7k 

We took into account that J = At'^/U. The magnitude of the t/U correction in ( pSj) is 
obvious, however one should be careful with the sign. To find it one needs to commute 
fermionic operators in order corresponding to Fig. 10b. The Dyson equation ( pT]) remains 
valid. So we can easily find the Green's function G^ , where index H indicates that the 
Hubbard model correction is taken into account. The plots of — - Im G^ {e, k) as functions 
of e for k = (7r/2, 7r/2), k = (7r/2, 0), k = (tt, 0), and k = (0, 0) are presented in Fig.ll. A 
plot of the quasiparticle residue Z^^ as a function of k along (1,1) direction is given in Fig. 
9. We see that the "Hubbard correction" causes the decrease of the residue outside of the 
magnetic Brillouin zone to be steeper. 



/CO / 7 \ "^ 7 

G'f(e,k)rfe^ l + -7k ~ 1 + 2-7k. (26) 

-oo \ t / t 



The sum rule (P3|) is changed. Now we have 

1. /•- ^u, ,. , f. J \^ , .J 

vr 

Let us comment on the definition (|l^) of the operator Cko-. Its normalization is adjusted for a 
system with strong antiferromagnetic correlations and it is close to that for rfko- (see eq.(^). 
However as a result the definition (0) differs from that usually accepted for a normal Fermi 
liquid by a factor y^- This is the reason why the sum rule ( ^6]) can be larger than unity. 
Generally the normalization can be chosen arbitrally. It is a question of convenience only. 
However, let us prove that the sum rule for the total number of electrons in the system is 



fulfilled. According to definition (|T6D 

cLck. = 2iVk,, (27) 



where 

(28) 




is the operator for the number of electrons. Due to the definition ( ]T7| ) of Green's function 
Gc one has the standard relation 

1 r°° 
- -Im / Gf,(e, k)de = (0|cLck.|0) = 2(0|iVk.|0). (29) 

TT J-oo 



Comparing with (^) we find 



|iVk.|0) = ^(l + 2^7k). (30) 



The operator for the total number of electrons is equal to 

N= J2 ^k.. (31) 

We put a "hat" to distinguish this operator from the number of sites N. Note that in all 
equations before we assumed summation over momenta inside the magnetic Brillouin zone. 
But in the eq. (^) we must sum over the full Brillouin zone. Finally from eqs. (PUD,(PTD 
one finds that the sum rule for the total number of electrons 

(0|iV|0)= E i(l + 2^7k)=iV (32) 

(T,kg/-u// 

is fulfilled. In conclusion of this discussion we would like to note that the origin of all these 
complications with normalization is very simple:The natural zone for the operator d^^ is the 
magnetic Brillouin zone. On the other hand the natural zone for Ck^ is the full Brillouin 
zone. This is the reason why one should be careful comparing these two operators. 

V. COMPARISON WITH EXPERIMENT 

Many of the experimental features observed are reproduced with the theory described 
here. The large dispersion between (0, 0) and (7r/2, n/2) and the assymetric quasi particle 
weight about (vr/2, vr/2) with the very strong dcrease in weight on moving beyond (7r/2, 7r/2). 
Also the lack of dispersion along (0, 0) and (tt, 0) as well as the very low quasiparticle weight 
is well reproduced. A major discrepancy between theory and experiment concerns the width 
of the quasi particle peak. The theoretical spectra (Figs. 8,11) have narrow peaks. On the 
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other hand the widths of the experimental ARPES spectraa are rather large, Fq ~ O.SeV ~ 
2.4J.Although the experiment has been done at higher temperature T = 350K this is 
probably too low to explain the widths of the peaks in terms of excited spin waves. It seems 
that another degree of freedom not considered in the present work must be of importance. 
One such possibility is the coupling to phonons. Such a coupling is expected to be quite 
large for the cuprates especially for the Cu-0 breathing mode since this strongly influences 
the stability of the Zang-Rice singlets. The coupling to phonons will result in a broadening 
of the quasi particle peak which can be described in terms of the Franck - Condon factors 
describing the probability that the system is left behind with a number of excited phonons 
upon the sudden removal of one electronEJ'Ej. Such a large width would indicate a strong 
coupling with phonons. In this paper we simulated the broadening with a Lorentz curve with 
a width Fq. The spectra obtained by convolution of — - Im (^^(e, k) with the broadening 
curve are given in Fig. 12. Agreement of these spectra with experimental onesEl is quite 
reasonable. 

VI. CONCLUSIONS 

In the present work we consider a modified t — J model at zero doping (insulating 
copper oxide plane) and zero temperature. We discuss the hole Green's functions with given 
pseudospin and given spin and derive a Dyson equation relating these two Green's functions. 
The Green's function with given pseudospin is very convenient for calculations, but it is an 
artificial object. In the real experiment one measures the hole Green's function with given 
spin. 

To describe the experiment we use hopping amplitudes for a modified t — J model as 
optained from LDA band calculations by Andersen et oB. Agreement of the theoretical 
spectra spectra with experimental onesa is quite reasonable after some broadening of the 
theoretical spectra.The physical origin of this broadening is not quite clear. Possible ex- 
planations are in finite temperature, and the contributions due to phonons and the electron 
phonon interactions. 
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Figure captions, 

Fig.l Plots of — - Im Gd(e, k) for different values of k. The Green's function is calculated 
in the self consistent Born approximation. 

Fig. 2 Quasiparticle residue Z^ of the hole Green's function Gd- Solid line corresponds 
to the direction to the Brillouin zone center: k G [(7r/2, 7r/2) — (0,0)]. Dashed line gives 
dependence along the face of the magnetic Brillouin zone: k G [(7r/2,7r/2) — (0,7r)]. 

Fig. 3 The vertices: a) - single hole creation, b) - hole + spin-wave creation, c) - usual 
hole-spin-wave vertex. Thick line correspond to Gc, and thin solid line corresponds to Gd- 
Dashed line is spin-wave. 

Fig. 4 Hole production mechanisms. Upper part of the figure describes wave function of 
initial Neel state: a - component without spin quantum fluctuations, b - component with 
spin quantum fluctuation. Lower part of the figure arises instantly from the upper one after 
removal of an electron with spin up. Part a does not contain spin flip, and it corresponds to 
the amplitude a^. Part b does contain spin flip, and it corresponds to the amplitude 6k,q- 

Fig. 5 Zero and first order in t diagrams for hole photoproduction. 

Fig. 6 Dyson equation relating Green's functions Gc (thick solid line) and Gd (thin solid 
line) . 

Fig. 7 An example of higher order correction which is taken into account in the Dyson 
equation (|2l|) . 

Fig. 8 Plots of — - Im G'c(e, k) for different values of k. The Green's function is found 



from Dyson equation (|2lD relating Gc with that in self consistent Born approximation. 



Fig. 9 Dependence of quasiparticle residues on k along (1,1) direction. 

Fig. 10 Mechanisms of an electron removal without spin flip: a) - t — J model, b) - Hub- 
bard model correction. 

Fig. 11 Plots of — - Im G^{e, k) for different values of k. The Green's function is found 
from Dyson equation (|2TD relating G^ with that in self consistent Born approximation. 

Fig. 12 Theoretical spectra for several values of k. Spectra obtained by convolution of 
— - Im, G^(e, k) with spreading imitated by the Lorentz curve. 
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